Nuclear fragmentation by tunneling 
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Fragmentation of nuclear system by tunneling is discussed 
in a molecular dynamics simulation coupled with imaginary 
time method. In this way we obtain informations on the frag- 
menting systems at low densities and temperatures. These 
conditions cannot be reached normally (i.e. above the barrier) 
in nucleus-nucleus or nucleon-nucleus collisions. The price to 
pay is the small probability of fragmentation by tunneling but 
we obtain observables which can be a clear signature of such 
phenomena. 

For a long time the problem of spontaneous fission has 
been of large interest for both experimentalists and the- 
orists Q . In fact spontaneous fission is a nice example of 
a quantum many body problem. Both the ./V-particle na- 
ture of the system, the atomic nucleus, and the quantum 
aspect, penetration of a barrier, has stimulated a huge 
literature production and still the problem is largely un- 
solved from a quantum mechanical point of view. To a 
minor extent, the "inverse process" of SF i.e. subbar- 
rier fusion has also stimulated many workers pl@]. It is 
the purpose of this letter to suggest a similar problem 
which involves the quantum iV-body features as well as 
the possibility of critical phenomena. 

Finite systems might show signatures which are typical 
of a phase transition. Experiments on fragmentation of 
atomic nuclei display typical features of a second order, 
power laws, critical exponents, etc.. or a first order phase 
transition (caloric curve). More recently similar features 
have been observed in metallic clusters and fullerens || . 
These are systems that in the infinite size limit have an 
equation of state that resembles a Van Der Waals. Finite 
size effects smooth divergences which become maxima for 
some function. Nevertheless careful analysis are able to 
extract to a good accuracy the values of critical expo- 
nents which are very close to the liquid-gas values. The 
problem of these experiments is that the role of dynam- 
ics is dominant [Q. As a consequence it is not possible 
to prepare the system at all desired excitation energies, 
densities or temperatures as for a macroscopic system. 
For instance in heavy ion collisions if the energy of the 
particle is too low we observe fusion-fission of the two 
nuclei, while, if too large, fragmentation occurs. Imag- 
ine that in some way we have been able to prepare the 
system at density p and temperature T. Because of the 
compression and/or thermal pressure, the system will ex- 
pand. If the excitation energy is too low the expansion 
will come to an halt and the system will shrink back. This 
is some kind of monopole oscillation. On the other hand 
if the excitation is very large it will quickly expand and 
reach a region where the system is unstable and many 



fragments are formed. It is clear that in the expansion 
process the initial temperature will also decrease. We 
could roughly describe this process with a collective co- 
ordinate R(t), the radius of the system at time t and its 
conjugate coordinate. Here we are simply assuming that 
the expansion is spherical. These coordinates are some- 
what the counterpart of the relative distance between 
fragments in the fission process Similarly to the fis- 
sion process we can imagine that connected to the col- 
lective variable R(t) there is a collective potential V(R). 
When the excitation energy is too low it means that we 
are below the maximum of the potential. That such a 
maximum exists comes, exactly as in the fission process, 
from the short range nature of the nucleon nucleon force 
(or the typical Lenard-Jones force for metallic clusters) 
and the long range nature of Coulomb. Thus, similarly to 
SF, we can imagine to reach fragmentation by tunneling 
through the collective potential V(R). When this hap- 
pens, fragments will be formed at very low T not reach- 
able otherwise than through the tunneling effect. The 
price to pay as in all the subbarrier phenomena is a low 
cross section. We call the fragments formed via tunneling 
Quantum Drops (QD). Of course a necessary conditions 
to make such fragments is to have enough energy and 
this could be easily estimated from the mass formula. 
For instance one could imagine to search for the events 
where 3 intermediate mass fragments (IMF) are formed 
(for which Z > 3) from an initial nucleus of mass A. 
^From these informations one can calculate the Q-value 
of the reaction and the corresponding excitation energy 
to obtain IMF. The closer the excitation energy to the Q- 
value the smaller the probability to observe the events. 
It is the purpose of this paper to show through a mi- 
croscopic simulation that this process is indeed possible 
and to discuss some features that undoubtedly charac- 
terize the formation of the Quantum Drops and that can 
be verified experimentally for instance in proton nucleus, 
gentle nucleus-nucleus or cluster-cluster collisions. 

The microscopic simulational study using Vlasov 
equation has been carried out for fusion and sponta- 
neous fission problems 0j below the Coulomb barrier. 
The dynamics has been formulated introducing a one- 
dimensional collective degree of freedom which develops 
in real- and imaginary-time. The method has reproduced 
well the fusion cross sections, fission kinetic energies of 
the fission fragments and other observables ||. The es- 
sential point is to introduce collective coordinates which 
enables us to treat the dynamics in imaginary time. 

Since the expansion of the system can be parametrized 
in a one-dimensional collective coordinate, we can, in 
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principle, apply the imaginary-time method to its study. 

In this paper we combine imaginary-time prescription 
with quantum molecular dynamics model to simulate 
fragmentation of finite system with relatively low exci- 
tation. 

InQMD @, each particle follows the equations of mo- 
tion given by the Hamiltonian. The effective interaction 
in the Hamiltonian used in this paper is taken from the 
QMD calculation of Ref. || which can reproduce the sat- 
uration properties of nuclear matter, binding energies of 
ground state nuclei and the energy dependence of the 
optical potential. To describe the sub-barrier dynamics 
of the system, we treat the collective radial expansion in 
real- and imaginary-time coupled to the motion of each 
particle. If we do not allow the tunneling, the system 
will shrink back at the turning point if there is one, or 
the system continues to expand if it has enough energy to 
overcome the attractive potential, i.e., there is no turning 
point. What we do here is to simulate the tunneling mo- 
tion when the collective variable reaches a turning point. 

We treat initially spherical finite system. Both the 
total center of mass coordinate and the total momentum 
are 0. The collective coordinate i? co ii and momentum 
P co ii for radial expansion is defined as 



-Rcoii — -7 ■ R», 

i 

i 

^ = Ri/\Ri\. 
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The equation of motion for each particle is coupled to 
the collective variables as follows. Before the collective 
radial motion reaches the turning point, each particle 
obeys the normal equation of motion, 



R 4 = dH/dPt, 
P, - -dH/dRi 



(4) 
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After reaching the turning point, the path of collective 
variables (real position and imaginary momentum) is 
generally obtained by the inverse potential method. We 
couple this inverse potential method for collective vari- 
ables and real-time evolution of each particle by intro- 
ducing the collective force F co a, 



R, = dH/dPi, 
P, = -dH/dRi 

F co n = Pcoll- 



2Pcoll&/A 
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Solving these equations of motion consistently, one ob- 
tains the real- and imaginary-time evolution of radial ex- 
pansion. 

When the collective expansion reaches the turning 
point again, we turn off the inverse potential and solve 
real-time equation of motion. If the collective variables 



never reach the turning point again, the event is not ac- 
cepted. Notice that if the system does not have enough 
energy to make fragments, then the second turning point 
will never be reached. 

To obtain any physical quantity, tunneling events have 
to be counted with their probability to happen, while 
events without tunneling can be counted as to have prob- 
ability 1. The corresponding probability II of tunneling 
event is calculated using the action S during the tunnel- 
ing process, 



n 



l + e 2S 



(9) 
(10) 



where i? a and Ry, are starting and ending points of the 
tunneling region, i.e., the two turning points. 
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FIG. 1. Snapshot of tunneling fragmentation. White cir- 
cles represent neutron and gray ones protons. The time in the 
simulation t is shown in each panel. After the system enters 
the tunneling region, the time is written as ti. 

We simulate the expansion of 230 U system. First we 
prepare the ground state of a nucleus and then com- 
press uniformly to an excitation energy E* from 5 to 
8 MeV/nucleon. Due to the fluctuations between events 
caused by the different initial configurations of the nuclei, 
the potential energy during the expansion is different for 
different events. Therefore the tunneling fragmentation 
occurs in some events where the potential energy is even- 
tually high, while there is no tunneling for events with 
lower potential energy. At sufficiently low E* , fragmen- 
tation occurs via tunneling only. In Fig. [I] we display 
snapshots of a typical tunneling event. The collective 
coordinate P(t) becomes zero at t = 100 fm/c. At this 
stage we turn to imaginary times as described above and 
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the tunneling begins. Notice that the system indeed ex- 
pands and its shape can be rather well approximated to 
a sphere at the beginning. But already at 300 frn/c (in 
imaginary time) due to the molecular dynamics nature 
of the simulation, the spherical approximation is broken. 
This shortcoming of our approach should be kept in mind 
because the calculated actions will be quite unrealistic. 
This is similar to the use of one collective coordinate (the 
relative distance between centers) in SF process. In that 
case many calculated features are quantitatively wrong 
but qualitatively acceptable. Since our Quantum Drops 
is a proposed novel mechanism we can only give qualita- 
tive features and the model assumption must be refined 
when experimental data will start to be available. 
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FIG. 2. (a) Reaction probability; (b) Average mass of 
large (Z > 3) fragments; (c) Average momenta of large frag- 
ments. Note that the horizontal axis is shown by the normal- 
ized excitation energy E* /Eo, where Eo = 7.0 MeV/nucleon 
(the threshold for M — 3 events.) 

In fact the calculated action S in the tunneling events 
for lower excitation energy is large and the probability 
of tunneling fragmentation is very small. Figure |^(a) 
shows the calculated probability of fragmentation with 



fragment multiplicity of 2 and 3 vs. E*/Eq. Here Eq 
is the threshold for having real events with M = 3. The 
value Eq depends on the detail of the force and in our 
case Eq = 7 MeV. This seems to us a quite unrealistic 
value; in fact when we compare our model to data on 
fragmentation (above the barrier) we underestimate the 
data, i.e. the model is less explosive. ^From the data, 
for instance Aladin data 0, we can estimate Eq « 3 
MeV/nucleon. 

For the imaginary events the probabilities arc obtained 
using the action S of the tunneling process. If only one 
type of reaction is considered the probability of both the 
real and the imaginary events are 1/2 at the threshold 
(barrier) energy. We see good agreement of probabilities 
for M = 2 near the vanishing point of real event. But the 
value of the probability is not 1/2 since there are many 
types of events we have to consider and the barrier height 
is fluctuating. For M = 3 events this effect is serious and 
the probabilities for real and imaginary events do not 
meet at any point. 
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FIG. 3. (Color) Dalitz plot for fragment mass. Red points 
indicate data from real fragmentation and blue are imaginary 
process. The excitation energy shown in each panel is nor- 
malized by threshold energy Eo- 

Although the probability of tunneling fragmentation is 
very small, one possible observable may be the fragment 
mass distribution. In Fig. ||(b) we plot the average frag- 
ment mass for intermediate mass fragment multiplicity 
M of 2 and 3. Although we see slight difference between 
imaginary and real events, the values are too close to dis- 
tinguish. Figure || is a Dalitz plot of three largest frag- 
ment mass in each event. Plotted with red color is frag- 
mentation without tunneling process and blue by tunncl- 
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ing. In this figure we can clearly distinguish both types 
of reactions. When the fragmentation occurs by tunnel- 
ing the size of fragments are more equally distributed 
(central region in Dalitz plot) than the case without tun- 
neling process. Equally distributed fragment size is not 
preferred by the Coulomb potential when position of frag- 
ments are concentrated in the early stage of fragmenta- 
tion. Such a situation is only possible with tunneling 
process. Without tunneling process, on the other hand, 
the low-energy fragmentation is only possible through se- 
quential decays. 

Another possible observable to distinguish the tunnel- 
ing events from the real events is the momentum distri- 
bution of produced fragments shown in Fig. ||(c). Again 
we see significant difference between imaginary and real 
events. It is rather natural that real events have larger 
value of fragment momenta since the eventual difference 
of kinetic energy of the expansion determines whether the 
system experience the imaginary process or not. How- 
ever, the difference of the process itself has some effects 
on the resulting fragment momenta: In the real events 
the fragments momenta are larger for cases with large 
fragment multiplicity M = 3 compared to M = 2. This 
is due to the different Q-value for events with M = 3 
and 2. In fact the difference of the total kinetic energy of 
large fragments for two types of events, which amounts 
to about 30 MeV, almost exactly corresponds the differ- 
ence of their Q-values estimated by a simple liquid-drop 
model. 

In imaginary events, on the contrary, fragment mo- 
menta for M = 3 and 2 are almost the same (within a 
numerical accuracy) though there are large fluctuations 
due to the poor statistics. This is due to the fact that 
below the threshold, the system must use all the exci- 
tation energy to form the final fragments. In order to 
do that the final fragments keep the lowest possible ki- 
netic energy. The final kinetic energy is essentially due 
to Coulomb repulsion at the second turning point. Thus 
precise experimental informations on the final kinetic en- 
ergies will give us an insight on the densities at freeze out. 
The comparison of fragment momenta between the frag- 
ment multiplicity M = 2 and 3 may give a signature of 
imaginary process in the fragmentation. 

In conclusion, we have simulated the fragmentation 
process by tunneling. The calculated probability of such 
phenomena is rather small, and the threshold energy for 
fragmentation in our calculation has been found to be un- 
rcalistically high. Therefore by employing more suitable 
interaction, the threshold energy and even the tunneling 
probability could change. Apart from the probability of 
the tunneling fragmentation, we have proposed two pos- 
sible observables for this phenomena. One is the Dalitz 
plot of produced large fragments. In tunneling fragmen- 
tation, fragment mass is more uniformly distributed than 
that of normal fragmentation. The other is the fragment 
momentum distribution, where we have observed that 
the fragment momentum are insensitive to the fragment 
multiplicity. 
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